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Abstract 

The comparative behaviour of normal and anomalous averages as functions of mo- 
mentum or energy, at different temperatures, is analysed for systems with Bose-Einstein 
condensate. Three qualitatively distinct temperature regions are revealed: The critical 
region, where the absolute value of the anomalous average, for the main energy range, 
is much smaller than the normal average. The region of intermediate temperatures, 
where the absolute values of the anomalous and normal averages are of the same order. 
And the region of low temperatures, where the absolute value of the anomalous aver- 
age, for practically all energies, becomes much larger than the normal average. This 
shows the importance of the anomalous averages for the intermediate and, especially, 
for low temperatures, where these anomalous averages cannot be neglected. 



Key words: Bose-Einstein condensation, dilute Bose gas, normal and anomalous aver 
ges, coherent matter waves 



PACS: 03.75.Hh, 05.30.Jp, 05.70. Ce 



1 



1 Introduction 



In the last years, the physics of dilute Bose gases has received much consideration both in 
theoretical and experimental aspects. An exhaustive list of references can be found in the 
recent review articles [1-4] and books [5,6]. The theoretical description of dilute systems 
with Bose-Einstein condensate is, to a large extent, based on the ideas of the Bogolubov 
theory, advanced in the papers [7,8] and thoroughly expounded in books [9,10]. This is 
because the present-day experiments deal in the majority of cases exactly with such dilute 
Bose gases. 

The starting idea of the Bogolubov theory is the breaking of gauge symmetry by means 
of the so-called Bogolubov shift for the field operator 

V>(r) = r](r) +^i(r) , 

where r/(r) is the condensate wave function, related to a coherent state [11,12] and defining 
the condensate density po(r) = |^( r )| 2 - As soon as the gauge symmetry is broken, in addition 
to the normal averages < ip^ (r)ip(r') > or < a\ak >, depending on whether the real or 
momentum space is considered, the anomalous averages < ip(r)ip(r') > or < a k a_ k > arise 
in the theory. 

The physical meaning and general features of the normal averages, which are directly 
connected with the reduced density matrices [13], seem to be familiar and clear to everyone. 
Contrary to this, the role and general behaviour of the anomalous averages are not always 
properly understood. There exists a very widespread delusion that the anomalous averages 
can be omitted at all, even at low temperatures. One often ascribes this unjustified trick 
to Popov, calling it the "Popov approximation". However, as is easy to infer from the 
original works by Popov [14-17], he has never done anything like that. He considered the 
temperatures in the close vicinity of the critical temperature T c . When the temperature 
T tends to T c , then the condensate density p tends to zero. The anomalous averages 
are proportional to p , hence, also tend to zero together with p . Contrary to this, the 
normal averages are proportional to the density of noncondensed particles pi, which is close 
to the total density p, when T ps T c . That is why in a narrow neighbourhood of T c the 
anomalous averages automatically become much smaller than the normal averages, without 
any additional assumptions. However, for a dilute gas at low temperatures T <C T c , the 
condensate density can become comparable with the total density, p ~ p, while, vice versa, 
the density of noncondensed particles can be much smaller than p. Then anomalous averages 
can be much larger than the normal ones, and there is no any grounds for omitting the former. 
Popov [14-17] has never neglected the anomalous averages at low temperatures T <C T c . 

It is the aim of the present paper to analyse the behaviour of the anomalous averages at 
different temperatures and to illustrate the properties of the anomalous averages as compared 
to those of the normal averages. Such a comparison reveals the role of the anomalous averages 
emphasizing the importance of their contribution, especially at low temperatures, where the 
anomalous averages can never be neglected. 

Throughout the paper, the system of units is employed, where the Planck constant h = 1 
and the Boltzmann constant ks = 1 are set to unity. 
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2 Normal and Anomalous Averages 



The behaviour of the averages can be rather different at varying temperature. It is possible 
to distinguish three principally different regions. For this purpose, the temperature T should 
be compared with the characteristic value p$ , where p = N/V is the particle density for N 
particles in volume V and 

$o = | $(r) dr (I) 

is the integral of the interaction potential 3>(r), assumed to be integrable. The value p$ 
characterizes the average strength of particle interactions. Another typical quantity is the 
critical temperature T c , which for a weakly interacting gas is close to the critical temperature 
for an ideal gas 
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where m is particle mass and £(3/2) = 2.612. The relation between p$o and T c can be 
obtained by introducing the scattering length a s through the equation 

$0 = 47T ^ . (3) 

m 

For the repulsive interactions, the scattering length is positive, which is assumed in what 
follows. Then we find 



p$ 



= 2 
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Keeping in mind a dilute gas, for which by definition 

pa 3 s « 1 , (4) 

we obtain 

(5) 

In this way, we can distinguish the following temperature regions. The interval of low 
temperatures is 

< T < p$ • (6) 
Intermediate temperatures correspond to the inequalities 

p$ < T < T c . (7) 

And the temperatures T ~ T c close to the critical point pertain to the critical temperature 
region. 

The critical region of temperatures in the close vicinity of the condensation temperature 
T c was analysed by Popov [14-16]. In this region, the anomalous averages turned out to be 
small as compared to the normal ones. This fact is very simple to understand, since when 
temperature tends to T c , the anomalous averages are proportional to the density of condensed 
particles p , which tends to zero as T — > T c . But the normal averages are proportional to 
the density of noncondensed particles pi, which tends, as T — > T c , to the finite value p of 
the total particle density. This is why in the close neighbourhood of T c the contribution of 
the anomalous averages is small as compared to that of the normal averages. 
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However at the temperatures outside the critical region, the situation can be drastically 
different. And we have to study these lower temperatures more attentively. We shall consider 
a uniform dilute gas, for which the inequality (4) is valid. 

In a dilute gas, for T <C T c , almost all particles are condensed. Then the Bogolubov 
theory [7-10] is applicable. Following this theory, one makes the Bogolubov shift of the field 
operators and retains quadratic fluctuations of the noncondensed particles, which results in 
the Hamiltonian 

H = -z iVoPo^o - piVo + u k a\a k + - A k (4 a -fc + O-feOfe) , (8) 

1 k^O 1 k^O 

in which N Q is the number of condensed particles, p = N /V is their density, 

k 2 

u k = — + p (<S> + <S> k ) - p , (9) 
2m 

fi is chemical potential, and 

A fc = p $fc , (10) 
with $fc being the Fourier transform of the interaction potential, 

$, = / $f r ) e - lkr dr . 



h = J $ (r) 



Under the considered conditions, N Q ss N is close to the total number of particles N. 

The Hamiltonian (8) is diagonalized by means of the Bogolubov canonical transformation 

a k = u k b k + v*_ k bl k , 

in which the coefficient functions can be chosen to be symmetric and real, such that u* k = 
U- k = u k and v k = v^ k = v k . They are defined by the relations 

2 2-, Afc 

u k - v k = 1 , u k v k = - — , 

2£fc 
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where the Bogolubov spectrum is 

£ k = y/ul-Al. (11) 
Under this transformation, Hamiltonian (8) is recast to the diagonal form 

H B = E + J2 £ kb{b k -fiN , (12) 



with the ground-state energy 



Eq = \ N oPo <5> + i J2(e k - uj k ) . 
1 1 k^O 
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The separation of the condensate contribution, as is well known, is meaningful only 
when the momentum distribution of particles becomes singular at the point k — 0, which is 
associated with the particle spectrum touching zero at k — 0. Only then there is a reason of 
separating the zero-momentum terms [18], which is equivalent to the Bogolubov shift of the 
field operators [9,10]. At the same time, this guarantees a gapless spectrum, in agreement 
with the theorems by Hugenholtz and Pines [19] and Bogolubov [9,10]. These requirements 
can be simply formulated as the condition 

lim e k = , e k > . (13) 
Then from spectrum (11), we get the chemical potential 

H = p % • (14) 
Another way of defining the chemical potential is by minimizing the thermodynamic potential 

VI = —T In Tr e~ pH 



with respect to A . Then from the equation 



an oh 

= < 77TT" > = , 
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dN 



in the frame of the Bogolubov approximation, one obtains the same chemical potential (14). 
Respectively, Eq. (11) acquires the known form of the Bogolubov spectrum 



£k — 



m 



2m, 



(15) 



The diagonal Hamiltonian (12) makes it possible to find explicit expressions for different 
averages [20-22]. Our concern here is the normal average 



n k = < a\a k > 



and the anomalous average 



(16) 
(17) 

whose absolute values are to be compared. 

First of all, we can derive a general relation between averages (16) and (17), following 
from the Bogolubov inequality 



o"fc = < a k a- k > , 



< AB > 



< (< AA + >< B + B >) , 



valid for the averages involving any two operators. Setting here A = a k and B = a^ k , we 
find 

| < a k a^ k > | 2 < (< a k a\ >< aL k a_ k >) . 

This yields 



o\ < n k (l+n k ) . 



(18) 
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Calculating expressions (16) and (17), we have the normal average 



n k = ^ + ^ coth (fh.) - I (19) 
2e k \2TJ 2 K ' 



and the anomalous average 

ITi- = — v , , , , , , 

2e k \2T 



a k = --^ coth ( % ) . (20) 



Inequality (18) holds true, since 

n k (l + n k ) -o\= 2sinh 

The spectrum (15) possesses the asymptotic properties 

e k ~ ck (k — > 0) , 

^ (fc -> oo) , (21) 
where c = ^ p § /m is the sound velocity, which tells us that e k varies in the range 

< e k < oo . 

It is, therefore, convenient to analyse the behaviour of Eqs. (19) and (20) with respect to 
the variable e k . 

At small excitation energy, such that 

e k < A fc , e k < T , 

the normal average (19) has the asymptotic behaviour 



TA k A k 1 I (^k_T___ A| \ £ 



nfc " 4 + 12T + 2A fe 2 + 1.3T A fe 90T 3 J 8A| ' (22) 
and the anomalous average (20) is 

TA fc A fc A k e\ 

ak ^~^f~ TzT + rm^ ■ (23) 

Hence at low energies or low momenta k, the values nfc and |<7fc| asymptotically coincide 
for all temperatures. The region of small momenta k usually gives the largest contribution 
to the integrals representing the observable quantities. Consequently, the anomalous average 
plays an important role. 

At large excitation energy, such that 

e k > Afc , e k > T , 

the normal average behaves as 

'Ht) 2 -(t) 4+ ^- (24) 



while the anomalous average, as 

* k ^-^(l + 2e-^\ . (25) 

This shows that for k — > oo, the absolute value of the anomalous average is always much 
larger than n k <C | cx fc | . 

From Eqs. (22) and (23), we have 

lim( W -n l ) = i(l- -|-) , (26) 

where we take into account that for the considered dilute gas po ~ p. Equality (26) demon- 
strates that, when T > p$o, then the value |<jfc| is smaller than n k at k = 0, though becomes 
larger than the latter for k — ■> oo, according to Eqs. (24) and (25). This defines the region of 
intermediate temperatures (7), where \a k \ is comparable with n k , being below n k for small /c, 
but surpassing n k for large A;. And for low temperatures from the interval (6), the absolute 
value of the anomalous average \a k \ is larger than n k in the whole range of the momenta 
k > 0. 

To illustrate pictorially the relative behaviour of the normal and anomalous averages, we 
calculate Eqs. (19) and (20) numerically. It is convenient to measure the temperature T in 
units of p$o an d to treat n k and a k as the functions of the dimensionless variable 

E = 4^ • (27) 



A 



k 



Then Eq. (19) can be rewritten as 



n(E) = ooth (I) - i ,28) 



<r(£) = - ^= coth f . (29) 



and Eq. (20) takes the form 

nl /••) 

2E V2T7 

The normal average (28) is non-negative for all E > 0, while the anomalous average (29) is 
always negative. To compare their absolute values, it is useful to introduce the notation 

A{E) = \a{E)\ (30) 

for the absolute value of the anomalous average (29). 

Figure 1 shows the typical behaviour of the anomalous value (30) and the normal average 
(28) for the low-temperature region (6), where A(E) is always larger than n(E). Figures 
2 and 3 demonstrate the changes in these functions, when the temperature rises to the 
intermediate region (7), where the averages are comparable with each other, so that A(E) < 
n(E) at small E, but at large E, this inequality inverts to n(E) < A(E). 

It is very instructive to picture the ratio n(E)/A(E), clearly showing how much different 
is the normal average n(E), as compared to the anomalous value A(E). Figure 4 illustrates 
this for the low-temperature region (6). As is evident, the anomalous average essentially 
prevails over the normal one. Figure 5 corresponds to the borderline between the regions of 
low and intermediate temperatures, where the anomalous average is yet dominant. Finally, 
Fig. 6 displays the ratio n(E)/A(E) for the intermediate-temperature region (7), when the 
inequality n(E) > A(E) for small E overturns to n(E) < A(E) for large E. 
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3 Conclusion 



In the theory with broken gauge symmetry, anomalous averages play an important role. Al- 
though in the critical region, close to the condensate temperature, they become much smaller 
by the absolute value, as compared to the normal averages, however at lower temperatures, 
the situation turns to be essentially different. At temperatures p<3> < T T c in the dilute 
Bose gas, the absolute value of the anomalous average is of the same order as the normal 
average. For very low temperatures T < p<3> , the anomalous average surpasses by its abso- 
lute value the normal average and can even be much larger than the latter. This analysis 
shows that, as soon as the gauge symmetry is broken, the anomalous averages come into 
play, being at low temperatures as important as the normal ones. It is, therefore, principally 
incorrect to omit the anomalous averages for T C T c . 

Note also the fundamental importance of the anomalous averages in the calculations of 
the particle fluctuations [23-25]. In the frame of the Bogolubov theory [7-10], the dispersion 
of the number-of-particle operator N is 

A 2 (A>) = jV[l + 21im(n fc + <7 fc )] . 

The long-wave limits of n k and a k separately are divergent, as follows from Eqs. (22) and 
(23). However, their singularities are the same in the absolute values while opposite in signs, 
because of which they cancel each other, yielding the finite limit 



limK + "*) = 2 ( 
As a result, the dispersion of N becomes 



mc 2 



TN 

A 2 (iV) 



mc z 



describing normal particle fluctuations. If one would omit o k in A 2 (A), this would lead to 
a divergent dispersion, hence, to a divergent isothermal compressibility k t = A 2 (A)/pTA, 
which would mean the system instability [23-25]. Thus, neglecting the anomalous aver- 
ages not merely can drastically distort quantitative values of observables, but also provokes 
qualitatively wrong unphysical consequences. 
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Figure 1: The absolute value of the anomalous average A(E) (solid line) and the normal 
average n(E) (dashed line) as functions of the dimensionless energy (27) for the reduced 
temperature T/p$ = 0.1 
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Figure 3: The same as in Fig. 2 for the reduced temperature T/p<3> = 10. 
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Figure 4: The ratio n(E)/A(E) of the normal average to the absolute value of the anomalous 
average as a function of the dimensionless energy (27) for the low temperature T/p$ = 0.01. 



14 



n(E)/A(E) 




2 4 6 8 10 



E 

Figure 5: The ratio n(E)/A(E) for the borderline temperature T/p$ = 1- 
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Figure 6: The ratio n(E)/A(E) for the intermediate temperature T/p<3> = 10. 
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Figure captions 



Fig. 1. The absolute value of the anomalous average A(E) (solid line) and the normal 
average n(E) (dashed line) as functions of the dimensionless energy (27) for the reduced 
temperature T/p$ = 0.1. 

Fig. 2. The anomalous value A(E) (solid line) and the normal average n(E) (dashed 
line) as functions of E for the reduced temperature T/p$ = 1.5. 

Fig. 3. The same as in Fig. 2 for the reduced temperature T/p$ — 10. 

Fig. 4. The ratio n(E)/A(E) of the normal average to the absolute value of the 
anomalous average as a function of the dimensionless energy (27) for the low temperature 
T/p$ = 0.01. 

Fig. 5. The ratio n(E)/A(E) for the borderline temperature T/p$ = 1. 
Fig. 6. The ratio n(E)/A(E) for the intermediate temperature T/p$ = 10. 
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